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ABSTRACT 


A  dclcrrninisiic  itiixiure  theory  is  presented  for  periodic  particulate  composites.  The  u'ode) 
is  constructed  by  introducin'!,  convenient  microdisplacement  and  microstress  variables,  and  by 
usiiiR  a  regular  asymptotic  lechniqut  with  multiple  stales.  Governing  equations  and  appropriatt 
boiu/dary  coiiditious  are  then  dt'duced  from  Heissner's  mixed  variational  principle  \2b  .  In  order 
to  test  the  atcurary  of  the  present  model,  harmonic  wavt-  propagation  is  examined  and  com¬ 
pared  with  a\ailal)le  experimental  data  for  CJIass /Epoxy  and  Steel/PMNl.A  composites  reported 
by  Kinra  and  Ker  16  .  Also,  the  effective  clastic  moduli  are  computed,  and  the  results  com¬ 
pared  with  other  analytical  metliods. 


►  •  .S.  .V 
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1.  INTRODUCTION 


Composiir  materials  exhibit  very  different  behavior  under  static  and  dynamic  loads,  as 
compared  to  homogeneous  isolro[)ir  solids.  Therefore,  proper  assessment  of  the  composite's  phy¬ 
sical  properties,  such  as  elastic  moduli  and  wave  dispersion  spectra,  is  important  for  most 
engineering  applicalions. 

Different  analytical  methods  for  the  determination  of  the  effective  elastic  moduli  of  particu¬ 
late  composites  have  been  compiled  by  Christensen  il'.  These  are;  (i)  the  variational  approach 
based  on  extremum  principles  |2.3.4,5.,  (ii)  the  probabilistic  approach  16.7,8.9  .  and  (iii)  the  self- 
consistent  scheme  approach  '10.11.12  .  At  this  moment  it  is  safe  to  conclude  that  formulas  are 
available  which  predict  accurately  the  effective  elastic  moduli  of  the  composites. 


In  order  to  clarify  dynamic  behaviors  of  particulate  composites  ultrasonic  measurements  of 
longitudinal  and  shear  wave  velocities  have  been  recorded  by  Kinra.  Petraitis  and  Datta  !l4  at 
low  frequency,  and  by  Kinra  and  .Anand  |15  at  long  and  short  wavelengths  for  random  particu¬ 
late  composites.  Later.  Kinra  and  Ker  jl6i.  again  using  ultrasonic  devices,  measured  longitudi¬ 
nal  wave  speeds  at  low  and  high  frequencies  for  periodic  particulate  composites.  Their  experi¬ 
mental  investigation  revealed  the  existence  of  stop  and  pass  bands  in  the  periodic  composites 
ll6  .  which  were  totally  absent  iti  the  random  composites  [14.15  .  To  date  there  is  no  contin¬ 
uum  model  which  accounts  for  such  wave  phenomena  in  particulate  composites  with  periodic 
microstructure. 

The  dynamic  behavit)r  of  random  particulate  composites  has  been  studied  theoretically  via 
harmonic  wave  propagation.  Mai  and  Bose  17i.  using  a  probabilistic  approach,  gave  values  for 
the  effective  wave  speeds  for  imperfectly  bonded  spheres  randomly  dispersed  in  the  matrix. 
Their  analysis  is  confined  to  the  long  wavelength  range.  Datta  ,18  obtained  longitudinal  wave 
velocities  through  tin  scattering  of  a  plane  P-wave  by  a  distribution  of  rigid  spheroids.  Later. 
Datta  !l9j.  using  a  probabilislit  approach  similar  to  that  in  T8  .  computed  longitudinal  and 
shear  wave  speeds  for  the  case  of  identical  elastic  ellipsoids.  His  results  19  arc  however 
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nondispcrsivc,  since  the  calculations  stern  from  a  long  wavelength  approximation.  Dispersive 
models  for  random  particulate  composites  have  been  proposed  by  Beltzer  et  at.  [20]  for  a  viscoe¬ 
lastic  matrix,  with  dilute  particulate  concentration,  and  by  Ben-Amoz  [21]  using  asymptotic 
expansions  to  derive  governing  equations  for  the  macro- motion.  The  latter  model  fails  to 
predict  the  experimental  data  reported  later  by  Kinra  and  Ker  |1C].  Gauriard  and  Ubcrall  [22] 
have  also  presented  a  dispersive  model  for  random  particulate  distribution  at  low  concentrations 
using  scattering  theory. 

A  deterministic  micromechanical  theory  is  proposed  here  for  the  simulation  of  harmonic 
wave  propagation  through  periodic  particulate  composites.  This  method  has  been  successfully 
applied  to  angle-ply  laminates  [23:  and  fiber-reinforced  composites  [24j.  It  consists  in  using  a 
regular  asymptotic  technique  with  multiple  scales  to  derive  appropriate  trial  displacement  and 
stress  fields.  This  procedure  yields  a  sequence  of  microboundary  value  problems  (MBV'P’s) 
defined  over  a  unit  cell,  which  characterizes  the  periodicity  of  the  composite  microstructure. 
Governing  equations  and  boundary  conditions  are  then  deduced  from  Reissner’s  mixed  varia¬ 
tional  principle  l25'.  The  lowest  order  model  of  this  MBV’P  method  is  equivalent  to  the  ’’0(1) 
homogenization  theory.”  by  Bensoussan,  Lions  and  Papanicolaou  [26]  and  Sanchez-Palentia  [27]. 
This  gives  the  effective  elastic  moduli.  However,  for  dynamic  problems  where  the  wave  length  is 
of  the  same  order  of  magnitude  as  the  inclusion’s  dimensions,  a  higher  order  model  is  necessary 
to  simulate  wave  dispersion  spectra.  Such  a  theory  is  named  as  0(t)  (at  least)  homogenization 
theory,  where  e  denotes  a  typical  ratio  of  micro-to-macrodimensions  of  the  composite.  To  date 
such  a  high  order  continuum  model  has  not  been  constructed  for  periodic  particulate  compo¬ 
sites. 

Propagation  of  harmonic  waves  through  periodic  particulate  composites  is  examined  using 
the  approach  outlined  above,  and  the  results  are  compared  with  the  experimental  data  recorded 
by  Kinra  and  Ker  [16]  for  Glass/Epoxy  and  StecI/^PNIMA  composites.  Also,  in  the  limit  when 
€  -•  0,  the  effective  elastic  moduli  are  computed  and  compared  with  Hashin’s  results  |2  .  In  all 


2.  FOKMULATION 


Consider  a  volume  V  whirli  contains  particulates  periodically  distributed  within  the 
matrix,  as  sliown  in  Fig.  1.  This  volume  is  referred  to  a  cartesian  coordinate  system  f,,  x-..,  ij. 

The  following  notation:  (  a  =  1.2  will  designate  quantities  associated  with  material  a. 
such  ifiat  n  ^  1  represents  the  partieuiate  and  o  -  2  the  matrix.  Barred  and  unbarred  quanti¬ 
ties  will  be  associated  with  ditnensiotial  and  notidimensiorial  variables,  respective!) .  Unless  oth¬ 
erwise  specified,  the  usual  cartesian  indicial  notation  is  employed  where  latiti  indices  range  from 
]  to  3  and  repeated  indices  imiily  the  summation  convention.  Also.  (  ),  and  (  )  ,  are  used  to 
denote  partial  differentiation  witli  respect  to  x,  and  time  T.  respectively. 

With  the  help  of  the  foregoing  notation,  the  governing  equations  for  the  displacement  vec¬ 
tor  v}°^  and  stress  tensor  associated  with  the  o*^^-cotistituent  are; 

a)  Equations  of  motion 


(1) 


whore  represents  the  mass  density: 

b)  Constitutive  equations  for  isotropic  constituents 


a-H  =  (2) 

where  A*“',  are  Lame's  constants,  the  infinitesimal  strain  tensor  and  6,^  the  Kronecker 
delta. 

c)  Strain-displacement  relations 


.  J.  L  l-'-  , 


(3) 


d)  Interface  continuity  conditions 


.tCI  --  sjf-l  -  /,  ^-1., 


on  A I 


where  .4/  denotes  the  pariiculale-mai rix  interface  surface. 


(4) 


»■)  Ihiiial  ronciitions  at  r  ()  and  a[)pr()f>rial(  boundary  rofidilions  along  the-  boundary  d  I  . 


Conditions  (a)  througlj  (f)  ronslitulf  an  initial  boundary  value  problem.  However,  due  to 
tli<  )arg«  nuinlxr  of  Ijcterogeneities  in  the  medium,  its  direct  solution  represents  a  formidable 
task.  The  purpose  ot  the  following  analysis  is  to  circumvent  this  difficulty  by  deriving  a  set  of 
p.iuial  differentia!  eijiiat ion^  with  constant  eoefficients.  whose  solution  will  thus  represent  an 
approximation  to  thc’  original  prolileni.  To  this  end.  the  basic  etjuations  are  nondimensionalized 
by  introduring  the  following  j>arameters: 

A  tyjiical  macrosignal  wavelength 

A  typical  particulate  sparing  or  cell  dimension 

C(„)  ,  reference  wave  velocity  and  macro  mass  density 

reference  modulus 

t|„)  =  A  typical  macrosignal  travel  time 

c  H  A  A  ratio  of  micro-to- macrodimension 

With  the  help  of  these  parameters,  nondimensional  variables  are  now  defined,  such  that 


£(r.) 


P[m) 


(5) 


Equation  (5)  shows  that  the  material  properties  arc  periodic  in  the  r, -space  in  which  the  perit)d- 
icity  can  be  characterized  by  the  cell.  This  representation  implies  that  stress  and  displacement 
fields  will  vary  according  to  two  ba.sic  length  scales;  (1)  a  macro  length  characteristic  of  the 
body  size,  and  (2)  a  micro  length  characteristic  of  cell  spatial  dimensions.  Further,  in  most 
applications  of  interest,  the  microscale  is  much  smaller  than  the  macroscale.  Therefore,  it 
apjitars  convenient  to  use  an  asymptotic  technique  based  on  the  two  scales  involved  in  the 
problem.  To  this  end.  new  independent  micro-coordinates  are  introduced,  such  that; 


•,’':i>'LV\VJVVtT.V 


As  n  consi'quciicc.  all  field  variables  now  brconie  functions  of  both  the  macrovariables  z,  and 
niicrovariabb's  z,  : 

F{z.  .  t)  ^  F  {i.  .  I,  ,  t  ;  <)  (7a) 

Sjialial  derivatives  of  sticb  a  fuiicliou  F  tlien  take  tb(-  form; 


,  J-i  .  <  ;  t)  -  -  ^{z,  .  ;  0 

(  ,» 


where  (  )  ,  is  used  to  denote  partial  differentiation  with  respect  to  x,  .  For  notational  sinipli- 
cit'  F  will  be  written  as  /'  in  what  follows. 

Applying  Eqs.  (5)  and  (7b)  to  Eqs.  (1)  -  (4)  yields  a  new  set  of  'Synthesized"  governing 
equations,  given  by: 


a)  Equations  of  motion 


b)  Constitutive  equations  for  isotropic  constituents 


<r}r'  =  -* 


c)  Strain-displacement  relations 


4  -  “i.t’  -  -  +  “'"11 

t  i  [  >.)  I 


d)  Interface  continuity  conditions 

^ j.V  on  A,  (11) 

The  periodicity  condition  implies  that  all  field  variables  will  as.sume  equal  values  on  opposite 
sides  of  the  cell  boundary.  Therefore  it  is  only  necessary  to  consider  a  single  cell  in  order  to 

determine  the  dependence  of  these  field  variables  on  the  microroordinales  x,  .  As  a  ron.sequence 
edge  laver  effects  are  not  included  in  the  model. 


Till*  const  met  ioti  of  tlic  present  mixture  model  is  accomplislied  by  resorting  to  Heissner's 
mixed  variational  principle  25  for  displacements  and  stresses.  This  variational  princiiile. 


applied  to  the  synthesized  fields  using  the  multiscale  representation,  takes  the  form; 


I 


h  f 


V'  l«' 


".'j’  -  '  -  ("'"V  -  “'“I 

(  i.j  ' 


-  el“>( 


dll  dr-.,  di: 


-  J  I  *7’.V  e,f'>)]c/.4 


diidr^dj^ 


1-1  *r 


V’  l<»1 


(12) 


-J. 


*  I'  i  ^ 

0  =  1 


\-  (o) 


dA 


^  I  ^  •  i  \  • 

where  V  denotes  the  i,  -domain  occupied  by  constituent  o  (see  Fig.  2),  a  is  used  for  the 

i' 

approximate  stress  field.  represents  the  traction  vector  on  the  surface  d  \'j  where  tractions 
are  specified  and  7"/  is  the  interface  traction  vector  defined  by 


r/=  on  A, 


(13) 


Also. 


^^1“’  --  7 


<<«,<“’  -I  -*  <5ul“U 

I  t  )■'  I 


(24) 


and 


1(0) 


2/1  (3A  -  2/y) 


rt  ^ 


-(o) 


(15) 


3.  ASYMPTOTIC  ANALYSIS 


The  asyiriptoiic  technique  starts  by  assuming  tlie  following  expansions  for  the  field  vari¬ 
ables: 


{u..r7  ,<;f )  -  S  .  o  J,. 7,  .f)  :  f  Cv  1 


(10) 


Substituting  (Hi)  into  the  governing  equations  (8)-(ll)  and  grouping  terms  in  equal  powers  of  i  . 
a  sequcrice  of  problems  is  obtained.  I'iie  first  of  these  yields 

(!■) 


«/<'*>,=  0  .  <7  ,  =  0 

ItW.I 


Equation  (17a)  implies  that  are  independent  of  tlie  microcoordinates  .  Together  with  the 
zerotli  order  expansion  of  (11a).  one  has  then 


“i^oj  ■  ^'i(0|(^id) 

The  remaining  sets  of  equations  obtained’ from  (8)-(ll)  are.  for  n,>  0: 


(18) 


_  (o)  _  -(o),,  (ol  y-  (o'  -  (o)  —  _ 

^  »  P  ‘'ilnl.K  '  ‘’}.(n).j  -  jii")  ‘’uln) 

|i(r»  1 


,  (o'  -  _  (o 

e,}(n)  -  ~  lu.'n 


■!"V;  ■*  “dn).  -  " 


I 


l).J  j(n  -■  )|.I 


(19) 

(20) 
(21) 


Equations  (19)-(21)  are  supplemented  with  interface  continuity  and  j  -periodicity  conditions. 
These  are.  for  n  >  0; 


and  e  are  r  -  periodic  on  cell  boundary 


(22) 

(23) 


Equations  (19)  through  (’2.3)  geiierate  a  series  of  niicroboundary  value  problems  (MB\  P's)  for 
and  .  which  are  to  be  solvi'd  sequentially.  Th<  first  of  the  MB\'P's.  corresponding  to 
r  ,'“!))  and  ti,)"); .  is  called  the  0(1)  MIU’P's  and  is  defined  by  (17b).  (19b).  (20).  (21).  (22b).  (23b) 
with  n  0  and  (22a).  (23a)  with  n  -  1. 
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Higher  order  terms  in  the  asymptotic  expansions  for  llie  displacement  and  stress  Helds' 
micros! run iires.  can  be  computed  with  suitable  integrability  and  normalization  conditions. 
However,  due  to  the  complexity  of  these  MBX’P's.  a  variational  apf)roach,  based  on  Reissner's 
mixed  principle  i25  is  preferred  here  in  order  to  construct  the  mixture  model.  This  construction 
can  be  carried  out  by  supplying  Eq.  (12)  with  appropriate  trial  displacement  and  stress  Helds 

which  must  satisfy  the  interface  continuity  condition  (22)  and  tlu'  i  -periodicity  condition  (23) 
on  the  cell  boundary.  .\s  a  result  only  approximate  solutions  of  the  MBN’P’s  are  necessary. 

These  trial  functions  are  obtained  by  modeling  the  cell  as  two  concentric  spheres,  as  shown 
in  Fig.  2.  where  [r,6,6]  are  spherical  coordinates,  such  that 

r,  =  r  sin^  cos9 
* 

Zj  =  r  sin^  sinS  (24) 

* 

Z4  =  r  cose- 

The  volume  of  the  cell  is  denoted  by  ,  its  boundary  is  defined  by  r  =  1  and  its  interface  sur¬ 
face  by  r  =  (n’b)'  *.  The  quantities  n*®*  denote  the  volume  fraction  of  material  q  and  satisfy  the 
relation: 

„(«)  ^  „(s)  -  j  ^25) 

In  terms  of  the  spherical  coordinates  (r.^.d).  the  z  -periodicity  condition  for  the  concentric- 
spheres  approximation  can  be  expressed  as 


F{r,6,d)  -  F{r.T  -f  ff.T  -f  ^) 


(20) 


4.  TRIAL  DISPLACEMENT  AND  STRESS  FIELD 

The  0(1)  stress  and  0(f )  displat  ement  fields  can  be  determined  by  solving  the  0(1) 
MHN'I’  s  wbicli  are  defined  by  (17b).  (lOb).  (20)-(2?.).  As  shown  in  |2C,.  these  MBN’P's  are 
excited  l)\  l,(o).,'  Therefore,  the  mixture  formulation  becomes  more  tractable  b)  introducing 
microdisplacement  variaiiles  wiiiclj  represent  t’,(o),,  -*  suO* 


where  V  is  the  volume  of  the  cell.  Since  is  excited  by  b’,|0)..  there  holds 


5.  =-  S,  (28) 

As  an  0(t )  trial  displacement  field,  the  following  form  tnay  be  used 


ti.W  -  5.  (»*,t)?,'“'(^i 


where 


(29a) 


t 


Considering  the  0(c*)  difference  of  the  average  of  u,*“^  over  V  Eqs.  (18)  and  (29)  yield  the 
following  trial  displacement  field 


«(“'(l*,I;  ,t.c)  =  l',f“>(l*.t)  +  ‘  U4^j(l*,xr.«) 


(30) 


where  is  given  by  (29).  Eqs.  (29)  and  (30)  show  that  the  mixture  displacement  variables 

are  t',"*,  l'}^\  S,  subject  to  the  constraint  (28).  In.scrting  (30)  into  (20)  and  (21)  with  n  =  0  and 
considering  the  O(i^)  differences  of  the  average  stresses,  the  0(1)  trial  stress  field  may  be  writ¬ 
ten  as 


(o) 

kj 


(31) 


As  with  displacement  variables,  an  0(c )  stress  variables  are  conveniently  introduced  according 
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to 

P.U^i)  ?  — f  O  lry'>dA  *  ^  -L-  f„  (32) 

t  V'  V'  *; 

>  ¥ 

h)U>gr«r>iig  l-Cj.  (8a)  over  \  and  iriaking  use  of  ihe  j  -periodkiVy  condition,  the  following 
mixture  etinalion  is  obtained 

P,  =  (33) 

where  the  average  operation  is  defined  by 


1 

riH  v  ' 


J  ,i)dj,  dz.,dii 

V  («l 


(34) 


From  (33)  it  is  seen  that  P,  plays  the  role  of  an  interaction  body  force  between  the  two  consti 
tuents  through  the  interface. 

As  an  0(< )  trial  stress  field,  the  following  form  may  be  used 


•')  -  /’»(**.*)  si“Vi  )<^., 

As  a  result,  the  trial  stress  field  may  now  be  written  in  the  form 


(35] 


k,«)  =  .0  +  <  ,t) 


(36) 


where  ci‘,%1  and  <7' )“(',]  arc  defined  by  (31)  and  (35)  respectively. 


f 


5.  MIXTURE  EQUATIONS 

Substituting  tin  trial  displficcmont  and  stress  fields  defined  by  (30)  and  (36),  respectively, 
iiito  Heissner’s  mixed  variational  j)rinciple  Eq.  (12)  and  using  Eqs.  (14)  and  (15).  the  following 
governing  equations  arc  obtained 
a)  liquations  of  motion 

"'“‘rjr.l  -  '  (-1)“'’^  =  (37a) 

-K  -  -4t  -  -F 

(*  '  n-  n '  ' 


where 


(>13,>1,)  =  E  yl<“>(iy“>) 


Also  arc  constants  that  depend  on  n*’*  and  n'**  only.  The  non-zero  values  of  these  con¬ 

stants  are  given  in  the  Appendix.  It  should  be  noted  that  in  Eq.  (3/),  the  superscript  a  used 
to  denote  "average"  has  been  dropped  from  rl,“'  and  r.'®*)  since  these  are  independent  of  the 
* 

microcoordinates  t,  . 
b)  Constitutive  equations 


’’Ul  9 


fW  -  A'”)  -  /'‘"’(f’..,  +  fV“') 

(39a) 

rirJ  -  A'“'  5.^,  ^  (S„  S.  +6,,  S,) 

(39b) 

r;(^)  _  rr(')  t 

P  ~  A,  '  '  ■+  A  f  Sn 

(39c) 

1  ) 

e..V.V 


•»_ 


r>V 

\*«r*j*  ■ 


.>r.'  r  ;  j".'  p,  v  j  m,  J  ,  i  }  Fj  1 1 1 1'!  >'1  ■  j  pj  ■  I  ,  ■  ■  ■  ij  i 
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^^  litrt' 


A>  B  ]/ 


} 

•>.  V 


/!'“> 


»  1  I  3A<"'  -,  -ill 


(») 


c)  Uiiundary  rondiiioriP  filoii;;  c)  V 


where 


SV}”'  -  0  or  ^  (  I)"'’’' 


iS,  -  0  or  AiPjti,  -  f. 


(40) 


(41a) 

{41b) 


*V/ 


. 

■K' 

/>;■ 

‘V 


A  aV 
/•*aV 


ri“^'=  _L  J  ri'*' 


»■  v^V. 


(42a) 


*•  •• 


T.^  —  i  f  ri“>  ?'">* *dz,*</T.* 

t>*  *1 


V  l.l 


(42b) 


v;  .■•■.' 


e)  Initial  conditions 


specify  .  5,  ,  S, ,  at  f  =  0 


(43) 


Equations  (37),  (39).  (41)  and  (43)  constitute  a  well-posed  boundary  value  problem  with  respect 
to  time  I  and  rnacrocoordinates  x,. 


•A  A  A  .* 


-  IG  - 


I::-;-:;: 


6.  EFFECTIVE  ELASTIC  MODULI 


The  effective  elastic  moduli,  which  relate  stress  averages  to  strain  averages,  are  dcTmed  b\ 
the  0(1)  homogenization  theory.  They  can  thus  be  obtained  by  letting  <  —  0  in  the  mixture 


Lq.  (37)  and  setting 


c.'-  =  r. 


Ktiuation  (37a)  v  ields 


a  I-'  =  U, 


where  ert"'.  v  (  l)-'r(;il 

0“  1  J 


s'"'  =  L  /)•“' 


while  Eq.  (37b)  yields 


'!!'  -  'S’  -  -iiT  'Si  -  -W  'Si  -  -  0 


This  last  equation  can  be  solved  for  the  S,  by  inserting  therein  Eqs.  (39a. b).  Substituting  in 


turn  these  expression.s  for  S,  into  (4Ga).  one  obtains 


<7  ('")  =  e*”’ 


where  ^  ~  o  o  jj,  o  jj,  o  ^  jj  ^12 


-  |b',  ,,  I  5  2.  t  jj,  (  5}  -■  {-'j2.  t  5.1  -  f  1,2-  {-12  -  {'Jl!' 


and  C'"''  is  the  effective  elastic  modulus  matrix.  Although  the  medium  is  heterogeneous,  it  is 
isotropic.  Consequently  only  two  itidependent  constants  should  suffice,  as  is  indeed  the  case. 


.S 


is  convcniciiilv  wriltcn  in  the  standard  form: 

A  2/i  A  A  0  0  0l>"> 

A  A  f  2/i  A  0  0  0 
A  A  A-t2;i000 

"  0  0  0  /i  0  0 

0  0  0  0  //  0 

0  0  0  0  0  f!  , 

Tlic  compositt  Lame's  constanis  A*”*  and  arc  given  in  liic  Appendix. 

In  order  to  assess  the  accuracy  of  llie  present  0(1)  mixture  theory,  the  results  for  tiic  elas¬ 
tic  constants  of  a  tungsten  carbide  cobalt  particulate  composite  are  presented  and  comjtared 
with  those  obtained  b\  Hashin  !2|.  Young's  modulus  E  and  Poisson's  ratio  i/  are. 

for  tungsteii  carbide  -  102  y  10*^  psi  ,  -  0.22  (51a) 

and  for  cobalt  fT*'*  =  30  x  10®  psi  .  r'***  =  0.30  (51b) 

The  variations  in  terms  of  the  particulate  volume  fraction  n’’’  of  the  composite  bulk  modulus 
and  shear  modulus  are  shown  in  Fig.  3,  where  very  good  agreement  is  observed  with 
Hashin’s  results  12;.  For  the  bulk  modulus  A'*”*,  the  two  curves  are  identical.  The  expre.ssion 
for  A'*"*!  derived  by  Hashin  |2]  represents  an  exact  result  when  the  medium  is  completely  filled 
out  with  concentric  spheres. 


7.  HARMONIC  WAVE  PROPAGATION 


In  order  to  lest  the  accuracy  of  the  present  0(t )  mixture  model,  phase  velocity  spectra 
have  been  compared  with  experimental  data  for  harmotiic  wave  propagation.  Due  to  the  com¬ 
posite  isotropn.  it  is  suffiriem  to  consider  harmonic  waves  propagatiiig  along,  say.  the  2,- 
direclioii.  Therefore  lei 


.  5.(2, ,tn 


{r.'"’  ,  e*’*'*" 


(52) 


where  f,'"  S,  rejiresent  constant  amplitudes.  Also,  k  and  denote  the  wave  number  and 
the  angular  frecptiency.  respteciivel) . 

Subsliiuling  (52)  into  (37)  and  using  the  constitutive  equations  (39).  yields  an  eigenvalue 
problem  for  (to.')  of  the  form: 

|ifl(/=  (tu)’:A/'  r  (53) 


where  L' =  ]  I/,'”  ,  t/r'  ,  5,' 


J  * 
S; 


s  * 
5,  , 


Cl"  ,  CJ" 


5, 


id) 


Cj=’ 


Sy 


r  *.T 
*5  1 


(54) 


\K'<  and  lA/  are  (12  x  12)  real  symmetric  matrices,  whose  elements  are  functions  of  the  volume 
fractions,  constituents  elastic  properties  and  wave  number.  Furthermore,  \M ,  is  a  diagonal 
matrix.  F'or  a  given  {>  k).  a  corresponding  eigenvalue  (tw)  can  be  computed.  The  phase  velocity 
Cp  is  then  determined  from 


(<*) 


(55) 


Iti  the  present  analysis  the  typical  cell  dimension  A  was  chosen  to  be  the  cell  radius  by  using  the 
concentric-spheres  model.  The  reference  elastic  modulus  and  density  used  for  the  nondirnetisioii- 
alization  arc: 


E 


(m) 


4 

o-  I 


^(m) 


r 

0=  1 


(5f>) 


The  dimensional  frequenev  /{fl-)  is  then  given  by 
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.\ui))rri(al  results  wtre  obtained  for  two  composites:  Glass/epoxv  and  Steel/PNfMA.  for  wliicb 
cxperiiiieiila!  data  have  been  recorded  by  Kinra  and  Ker  jl6  for  longitudinal  waves.  The  coti- 
stituents'  physical  properties  are  given  in  Table  I.  in  which  the  values  for  Poisson's  ratio  are 
estimates,  lii  the  simulation  A  was  computed  from  the  particulate's  radius  which  was  I  rnti  for 
glass  and  0.5'>  mm  for  steel. 

Phase  velocity  profiles  for  each  composite  are  shown  in  Figs.  4  and  a.  where  two  distinct 
modes  are  displayed;  a  longitudinal  (P-  mode)  wave  and  a  shear  (S-  mode)  wave. 

The  first  two  branches  in  the  spectrum  for  the  P-  mode,  which  are  believed  to  be  suffi¬ 
ciently  accurate  within  the  present  calculation,  are  plotted  as  functions  of  frequency.  Figures  4 
and  5  clearly  bring  out  the  dispersion  characteristic  of  the  composite.  The  experimental  data 
correspond  to  the  longitudinal  mode,  and  it  can  bo  seen  that  the  correlation  with  the  present 
results  is  quite  satisfactory.  In  their  paper.  Kinra  and  Ker  !l6,  have  predicted,  at  higher  fre¬ 
quencies  where  the  wavelength  approaches  the  particulate's  radius,  additional  second 
(f=1.4MHz)  and  third  (f=-2.(J  Mllz)  stop  bands  for  Glass/Epoxy,  and  additional  third 
(f=].5  MHz)  and  fourth  (f=1.9  MHz)  stop  bands  for  Steel/PMM.A.  Their  experiments  however 
did  not  reveal  these  higher  stop  bands.  They  have  conjectured  'Yhai  this  is  probably  due  to  the 
excitation  of  re.sonances  of  the  spherical  inclusions  which,  at  higher  frequencies,  tend  to  dom¬ 
inate  the  resonances  of  the  unit  cell."  Figure  4  (Glass/Epoxy)  shows  an  additional  stop  band  at 
f=1.53  MHz.  which  is  in  fairly  close  agreement  with  the  value  predicted  by  Kinra  and  Ker  ;1G  . 

Finally,  for  completeness,  the  first  branch  of  the  shear  wave  velocity  spectrum  is  also 
shown  in  Figs.  4  and  5.  for  each  composite.  Due  to  the  composite  isotropy,  the  vertically  (.SV ) 
and  horizontally  (SH)  polarized  wave  speeds  are  identical. 


8.  CONCLUSION 


A  microriicclianical  niixiurp  model  for  i>oriodic  paniculalc  composiies  was  presented  based 
on  a  two-scale  asymptotic  expansions.  Governing  e‘C|Uations  and  appropriate  boundary  condi¬ 
tions  were  deduced  from  lUissner's  mixed  variational  principle  .25  .  Expressions  for  the  effective 
elastic  moduli  were  derived  and  compared  with  llashin's  results  2  .  Also,  haririonic  wave 
dis{)ersion  sjxxtra  were  obtained,  and  a  good  correlation  with  the  experimental  data  of  Kinra 
and  Her  IC  for  Glass  Epoxy  and  Steel  I’MMA  composites  was  observed. 
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